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Quantum walks, both discrete (coined) and continuous time, form the basis of several recent 
quantum algorithms. Here we use numerical simulations to study the properties of discrete, coined 
quantum walks. We investigate the variation in the entanglement between the coin and the position 
of the particle by calculating the entropy of the reduced density matrix of the coin. We consider 
both dynamical evolution and asymptotic limits for coins of dimensions from two to eight on regular 
graphs. For low coin dimensions, quantum walks which spread faster (as measured by the mean 
square deviation of their distribution from uniform) also exhibit faster convergence towards the 
asymptotic value of the entanglement between the coin and particle's position. For high dimensional 
coins, the DFT coin operator is more efficient at spreading than the Grover coin. We study the 
entanglement of the coin on regular finite graphs such as cycles, and also show that on complete 
bipartite graphs, a quantum walk with a Grover coin is always periodic with period four. We 
generalise the "glued trees" graph used by Childs et al. [STOC, 59, (2003)] to higher branching rate 
(fan out) and verify that the scaling with branching rate and with tree depth is polynomial. 
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I. INTRODUCTION 

One of the most important tasks on the theoretical side of quantum computing is the creation and understanding 
of quantum algorithms. The recent presentation of several quantum algorithms based on quantum versions of random 
walks is particularly important in this respect, since they provide a new type of algorithm which can show an 
exponential speed up over classical algorithms, to add to those based on the quantum Fourier transform. Childs et 
al. [l| have produced a scheme for a continuous time quantum walk that can find its way across a particular "glued 
trees" graph exponentially faster than any classical algorithm, while Shenvi et al. [3] proved that a discrete quantum 
walk can reproduce the quadratically faster search times found with Grover's algorithm for findinga marked item 
in an unsorted database. Generalizations to finding subsets of items have also been developed [s',^, providing 
polynomial speed up over classical algorithms. For an overview of the development of quantum walks for quantum 
computing, see the recent reviews by Kempe Q and Ambainis These results are extremely promising, but still 
a long way from the diversity of problems for which classical random walks provide the best known solutions, such 
as approximating the permanent of a matrix [8;], finding satisfying assignments to Boolean expressions (/cSAT with 
/c > 2) (9|, estimating the volume of a convex body and graph connectivity Classical random walks underpin 
many standard methods in computational physics, such as Monte Carlo simulations, further motivating the study of 
quantum walk algorithms. 

Like classical random walks, quantum walks come in both discrete time [l2l . [isl [l^ . [isj . and continuous time (l6j 
versions. The discrete and continuous time versions of classical random walks can be related in a straightforward 
manner by taking the limit of the discrete walk as the size of the time step goes to zero. In the quantum case, the 
discrete and continuous time walks have different sized Hilbert spaces so there is no simple limit that relates the two 
basic formulations. There is also an example of a problem where the algorithmic powers of discrete and continuous 
time walks differ. Spatial search, where there is a cost associated with moving from one data element to another, can 
be accomplished faster with a discrete time quantum walk, but a continuous time quantum walk only performs as 
well for spatial dimensions greater than four [l7|. A continuous time walk with extra degrees of freedom has also be 
formulated by Childs and Goldstone [18| that does correspond to the limit of the discrete time walk and can perform 
equally well on spatial search. 

Our work in this paper investigates the pr oper ties of coins in discrete quantum walks. We follow on from prior 
work on quantum coins bv lMackav et al.l 19j and iTregenna et "al] 20| . broadening the types of graphs studied. The 
question of what is particularly quantum in a quantum walk is an interesting one which has attracted much attention 
[2l|,[22,[2l]. In this paper we address this issue by investigating the evolution of the quantum mechanical entanglement 
as the quantum walk progresses. We quantify the entanglement between the coin and position for example in a coined 
walk by using the von Neumann entropy and show how the entanglement oscillates and approaches asymptotic values 
depending on the choice of initial state and coin bias. 

The paper is organised as follows: Walks on infinite lattices are discussed first, starting with the simple walk on a 
line in Sec.[Tll and progressing to walks on lattices in two spatial dimensions in Sec. IIIII Section HVl considers walks on 
finite graphs, including the iV-cycle, cycles with diagonals and complete bipartite graphs. We then consider quantum 
walks on the "glued trees" graph of and generalise it to higher branching rates in Sec.|Vl Finally, we summarise 
and conclude in Sec. |VT1 

II. WALK ON AN INFINITE LINE 

In a classical random walk on a line, a particle moves either left or right according to the state of a classical coin 
where heads means right and tails means left (or vice versa). For a quantum version of a random walk, the coin is 
a qubit that can be in a superposition of heads and tails, so the particle moves left and right into a superposition of 
positions. This evolution of the walk is governed by a coin operator that acts on the quantum coin at each step of 
the walk. 



^(gcn) _ / ^ \/l - Pe"^ \ N 
^2 - ^Tl—pe^^ -^e''(«+<^) ) ' 

where < 6*, < tt are arbitrary angles, < p < 1, and we have removed an irrelevant global phase so as to leave the 
leading diagonal element real. Equation ([1]) represents the most general expression |24j ] for a unitary coin operator 
with two degrees of freedom. In this expression, the factors p, and 6, cf) determine the bias and the phase angles of 
the coin respectively. If we set p = 0.5 and 6 = cj) ~ the following expression, called the Hadamard coin operator is 
obtained: 
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FIG. 1: Probability distribution for a quantum walk on the line after 5000 steps, using a Hadamard coin, equation ([2]), and 
symmetric initial state (j-R, 0) + ijL, 0))/\/2- Only even positions shown since odd positions are unoccupied. 




This is an unbiased coin operator, as it chooses the directions left and right on a line with the same probability. We 
label as \L) and \R) the basis states of the coin, which can correspond to spin-up and spin-down states respectively. 
We denote the position on the line by so the joint state of a particle at position x with a coin in state \L) can 
be written \L,x). For the quantum walk on a line, the phases in the general coin operator {9 and 4>) appear in the 
evolution of the walk only in the combination {9 + tp), so as shown by Bach et al [2J|, their effect is equivalent to 
varying the phase (3 in the initial coin state |?Ao)- 



IV'o) = V^|i> + v/T^e"3|i?), 



(3) 



where "q is the bias in the initial state, and [3 is the relative phase between the two components. This leaves only the 
bias p in the coin operator affecting the outcome of the quantum walk on a line, and, without loss of generality, we 
can consider coin operators of the form 



.i(bias) 



\/i - p -Vp 



(4) 



After "flipping" the coin with the coin operator, the particle moves to adjacent positions according the the coin state; 
this is expressed mathematically as a conditional shift operator 



S\L,x) = \L,x-l) 
S\R,x) = \R,x + l) 



(5) 



One complete step of the quantum walk is thus given by the unitary operator U = S{C (8) 1). 

The position probability distribution of a quantum walk on a line is by now well-known, an example with a 
Hadamard coin operator, and initial state of -^(1^) + ^ |0)^ after 5000 time steps is shown in figure [TJ 



A. Entanglement between coin and position 



Since the quantum walk dynamics are unitary, the system remains in a pure state and we can use the entropy of 
the reduced density matrix of the coin to quantify the entanglement between the coin and the particle's position, 

Ecit)^-J2x,\og^{X,) (6) 
i 

where {A^} are the eigenvalues of the reduced density matrix of the coin at time t (in the case of the walk on a line 
there are just two eigenvalues). Figure [5] shows how Ec{t) varies for different initial coin states using a coin operator 
equation (|4]), with bias p = 0.2. This shows that the entanglement approaches a limiting value that varies between 
zero and one depending on the initial state of the coin. The rate of convergence to the limiting value also depends 
on the initial state of the coin, with symmetric initial states (77 = 0.5) converging fastest (i.e. oscillations about the 
asymptotic value die away fastest). 
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FIG. 2: (Left) Entanglement Ec{t) for a walk on a line with a biased coin, equation (UJ with p — 0.2, and two different initial 
states, asymmetric \R,0) (blue) and symmetric {\R,0) + i\L,0)) /\/^ (green) [23]. (Right) Entanglement Ec{t) for a walk on a 
line with coin operator bias p = 0.5 (a Hadamard coin), and the same two initial states. 
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FIG. 3: Entanglement Ec{t) for a walk on a line: asymmetric initial coin state \L) shown left, and symmetric initial coin state, 
(|L) + i\R))/V^ shown right, for varying coin bias p, over 20 time steps. 



B. Limiting value of the entanglement 

For the unbiased (Hadamard) coin operator^ (p = 0.5), whatever initial coin state \ipo) is chosen, the asymptotic 
value of the entanglement Ec{t — > 00) — >■ £Had — 0.872. Note added: this statement is not correct, the asymptotic 
value varies with the initial coin state (26l . [IJ. However, the rate of convergence is very different for different initial 
coin states, with more symmetric initial states converging faster, see figure [2] 

For biased coin operators, the picture is more complicated. We have studied the limiting value of the entanglement 
for two different initial state, \L, 0) (asymmetric) and {\L, 0)+i\R, 0))/\/2 (symmetric), see figure[3] In the asymmetric 
case the entanglement converges to a limiting value for all p. The limiting value of the entanglement Eooip) increases 
monotonically from (p — 1) to 1^ {p — 0^) and is discontinuous at p = 0. For p — the coin operator becomes 
the Pauli spin operator ax and the entanglement is zero for all time steps, showing the coin and particle remain 
disentangled. In the symmetric case the entanglement converges to a limiting value for all except the extreme case of 
p = where the entanglement oscillates between the minimum and maximum values and 1, as can easily be verified 
analytically. The asymptotic value of the entanglement Eoo{p) increases to 1 as p increases. The other variable factor 
is the period of the oscillations about the convergent value. In both cases, except for p = this period increases to 



Roldan, Knight and Sipe have also studied this case analytically (unpublished). 
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FIG. 4: Amplitude of oscillations of entanglement versus initial coin state parameters a and fi, for a Hadamard coin operator 
applied for 200 time steps. 

infinity as p ^ 1. As noted for the unbiased (Hadamard) coin, the rate of convergence to the asymptotic entanglement 
is faster for the symmetric initial coin state. 

C. Rate of convergence 

In order to quantify the rate of convergence of the entanglement to its limiting value, we considered the magnitude 
of the entanglement at a fixed time while varying the initial state. It is convenient to write the initial coin state as 

IV^o) = cos(a)|i?) + sin(a)e*'^|L), (7) 

where 77 = cos^(q;). Figure [4] shows how the amplitude varies with both a and /3. First consider the case where a = ^ 
for r\ = 0.5 i. e. a symmetric initial state. We find that the entanglement is proportional to cos^(/3) along the dotted 
line in figure [H We can also fit a formula for the minimum, 

a = Jcos(/3) + J. (8) 

This is the white line in the blue region in figure [4] where the entanglement oscillations are smallest, i.e., fastest 
convergence. 

To summarise the results for the quantum walk on a line, we find the various behaviours of the entanglement are 
governed as follows. The asymptotic value E^a reached by the entanglement is a function of both the coin bias p and 
the initial state |'0o)- For ^ fixed number of time steps t, the period of oscillation of the entanglement Ec{t) around 
Eoo is a function of p only. For the special case of p = 0.5, an unbiased coin, E^o = E^ad has the same value of 
0.872 . . . for all choices of initial coin state |V'o)- 

III. LATTICES IN TWO SPATIAL DIMENSIONS 

A. Higher dimensional coins 

For lattices with more than two edges meeting at each vertex, there is a far wider range of unitary coin operators, 
since the coin must now have as many degrees of freedom as there are choices of path. Since the range of higher 
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dimensional coin operators is too large for systematic numerical stu dy, for the remainde r of this paper we concentrate 
on two natural choices. The Grover operator was first introduced bv lMoore and Russelll 28] in their study of quantum 

walks on the hypercube. Based on Grover's diffusion operator, it has elements [C^^^]i,j = | — Sij, i.e., 



.(G) 



(2/d-l 21 d 
21 d 2/d-l 

\ 2/d 2/d 



For example, the d = 3 case is 



2/d \ 
2/d 



2/d - f / 



(9) 




(10) 



Except in the c? = 4 case, the Grover coin is biased, since the incoming direction (corresponding to the diagonal entry) 
is treated differently from the outgoing directions. However, it is symmetric under interchange of any outgoing coin 
directions, and is in fact the symmetric unitary operator furthest from the identity. The d ~ A Grover coin is the only 
unbiased Grover coin since all the entries are ±i 



.(G) 




(11) 



The DFT (discrete Fourier transform) coin is unbiased for all d, but asymmetric in that you cannot interchange the 
labels on the directions without changing the coin operator: each direction acquires its own phase shift. For d = 3, it 
looks like 



1 / 1 1 1 



(12) 



where = e^*'^/^ and w| — e^^iTr/s ^^j-g -j-j^g complex cube roots of unity. The d-dimensional DFT coin can be written 
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(13) 
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where ld is the complex d'th root of unity. 



B. Cartesian lattice 



For a 2-dimensional Cartesian grid, there are four edges meeting at each lattice site, so a d = 4 dimensional coin is 
required. The quantum walk is a generalisation of the walk on a line. We tested both Grover and DFT coins (d = 4 
versions) and found a similar range of behaviours for the entanglement between the coin and the position as for the 
walk on a line, only compounded by having twice as many directions. So, for example, the period of the oscillations 
about the asymptotic value is now a more complicated pattern of two frequencies. We looked for a correlation between 
the rate of convergence of the entanglement and the degree to which the quantum walk spreads out over the lattice. 
Most choices of initial state for the Grover and DFT coin operators produce a high probability of finding the particle 
on or near the starting point, with only one special initial state giving a high rate of spreading, compare the two 
distributions in figure [HI taken from |20|. Spreading is a property of random walks that can be useful for efficient, 
uniform sampling, compare [29j . The entanglement converges much faster for the quantum walk that spreads out 
in the ring, see figure [6l The entanglement between the coin and the position thus provides a way to monitor the 
progress and character of the walk. 



FIG. 5: Typical (left) and ring-like (right) spreading distribution for a Grover coin on a 2-dimensional Cartesian lattice (from 
[2^). Axes represent position {x and y) with the z axes indicating the probability of finding the particle at that position. 
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FIG. 6: Entanglement Ec{t) for the spike (blue) and ring (green) distributions on a 2D Cartesian lattice, using a Grover coin. 
A DFT coin operator produces a very similar result. 



C. Triangular lattices 

Higher dimensional lattices that lie in a plane (two spatial dimensions) can be constructed in a number of ways. 
We studied two examples: a tessellation of equilateral triangles produces a lattice with d — 6; and adding diagonals 
to a Cartesian grid, makes a "first and second nearest neighbours" lattice with d ~ 8. These are illustrated in figure 

m 

We tested a number of different initial states with both Grover and DFT coins, again looking at correlations between 
the amount that a walk starting from the centre of the grid spreads and the oscillations in the entanglement. To 
quantify the spread we studied the mean square deviation of the probability distribution from the uniform probability 
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FIG. 8: Probability distribution after 49 steps of a quantum walk on a d = 6 grid using a Grover coin operator: (Left) starting 
with coin state X^iLi I*)' point (50,50), and (Right) starting with coin state X^f=i 1^) ^ I*)- Axes 

represent position in the a;j/-plane and probability as in figure [S] 



distribution 

a^^'w = E(|'/'(-'*)I'-]^)'' (14) 

where r is a lattice site in the set it is possible to reach after t steps, and N{t) is the number of such lattice sites (so 
1/N{t) is the average probabihty per site). If the walk is spread out evenly over the lattice then AP (t) will be small, 

whereas if it is concentrated on parts of the lattice, AP (t) will be larger. 

For the both the d — 6 and d — 8 grids, as on the rectangular grid, the Grover coin can produce two kinds of 
behaviour: fast spreading distributions and distributions concentrated nearly all close to the origin, depending on the 
choice of initial states. The amplitudes of the oscillations in the entanglement decrease quickly for fast spreading but 
only slowly for the distributions stuck near the starting point. This is illustrates in figures [5][ni 

In contrast to the Grover coin operator, the DFT coin produces good spreading for almost all chosen initial states 
and the entanglement converges faster for these cases. To investigate the correlation between the entanglement and 
the spreading of the walk more thoroughly, we wrote the initial coin state as 

For d = 6, we fixed one of the phases and varied the other d ~ I over all the permutations, allowing repetitions of 
(j>j — Trj/d, {j — 1,. . . ,d, total of 6^ values, thus) to see the interference effects on the spread. For the Grover coin 

operator, most of these states give a large AP^(i) while for the DFT most give a small AP^(t): the average (over 

these initial states) after t = 15 steps is (AP^(t)) = 0.0501 for the Grover and (AP^(t)) = 0.0015 for the DFT walk. 



10" 



10" 



.-4 



10 20 30 40 50 
Time step 



1.4 
1.2 
1 



10 



20 30 40 50 
Time step 



FIG. 9: Spreading as indicated by the mean square deviation from the uniform distribution AP (t) (left) and entanglement 
Ec{t) (right) for the cases shown in figure |8] spike (blue) and ring (green). 



For the minimum AP (t) in this set, the entanglement oscillations decay fastest, as seen in figures |9] (for d = 6) and 

[TT] (for d = 8). The minimum values of AP (t) occur for the initial states = (0,0,0,0,0,0) for the Grover coin 
operator and = (0, 0, 4, 1, 1, 5) for the DFT coin operator. 

For d = 8, we fixed one of the phases and varied the other d — 1 over only the following subsets of all the arrays 
of 0j = TTj/d: (i) j e 0,d- 1, = (0,0,0,0,0,0,0,0),(0,0,0,0,0,0,0,d- 1), . . . , (d - 1, d - 1, . . . , d - 1, d - 1, 0); 
(n) permutations of (0, 1, 2, . . . , d — 1) (without repetitions) to see the interference effects on the spread. Although we 
didn't vary over as large a range as in the d = 6 case, we could observe a similar behaviour: better spread for DFT 
against Grover, entropy converging quickly for better spread. 

The average (over these initial states) after t = 15 steps is (AP^(t)) = 0.0508 for the Grover and (AP^(t)) = 0.0027 
for the DFT walk. (For the d — 8 grid N{t) = {2t + 1)^, since there are diagonals everywhere). The minimum 

values of 7V(i)AP (t) occur for the initial states — (0,0,0,0,0,0,0,0) for the Grover coin operator and {(pj} — 
(0,7,7,7,7,0,0,7) for the DFT coin operator. For the d = 8 grid we also tested Hadamard coins (a simple 2D 
Hadamard coin for each of the 4 pairs of opposite directions) and found they displayed similar properties to the DFT 
coin. 

Thus, for all tested coin operators and coin initial states, we have verified that those walks with good spreading also 
show small amplitudes of oscillation in the entanglement of the coin, pointing to a quick convergence. In particular, 
the DFT and Hadamard coin operators produce faster spreading than the Grover coin for most choices of initial state 
on these lattices of higher degree. This can be explained by noticing that for d > 4, the Grover coin is biased so that 
it favours returning along the edge that it arrived from. This will tend to reduce its spreading power. The DFT coin 
is unbiased, so its spreading power is affected only by how the different phases cause interference effects. For larger d, 
there are more different phases, so less opportunities for them to all cancel out. The Hadamard coin is also unbiased, 
and it does not mix between the different orientations of the pairs of edges, so on these triangular lattices it produces 
a spreading equivalent to the spreading on a line. 



IV. WALKS ON FINITE REGULAR GRAPHS 



We now turn to quantum walks on graphs with a fixed number of vertices so the walk is bounded and the notion of 
spreading is no longer the relevant property. Quantum walks on finite graphs were first investigated by Aharonov et 
al [14| . who showed that while the instantaneous distribution of a quantum walk on these graphs does not converge 
(being unitary and reversible), a suitably defined time-averaged distribution always converges, though this distribution 
need not be uniform (in the classical case the limiting distribution is always uniform). The interesting questions are 
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FIG. 10: Grover coin on d = 8 grid, after 39 time steps, start with coin state X^f^i 1*) ^ !*)• ^.xes represent 

position in the xy-plane and probabihty as in figure [5] 
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FIG. 11: Spreading as indicating by the mean square deviation from the uniform distribution AP (t) (left) and entanglement 
Ec{t) (right) for the DFT coin operator on d = 8 grid, after 39 time steps, starting with coin states Yli=i K) (blue) and 

^Ef£!N)+i;^Et./.+iN) (green). 

thus how fast the quantum walk converges to the Umiting time-averaged distribution, and whether the distribution is 
uniform. We are also interested in whether the walk shows periodic behaviour in the instantaneous distribution [20j . 
and if so, under what conditions. 



A. Ai'-cycles 



The quantum walk on a line can be converted to an A'^-cycle by taking a line segment of length N , and applying 
periodic boundary conditions. Clearly, for t greater than N/2, when the walk starts to wrap around on itself, the 
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FIG. 13: Comparison of time-averaged entropy (thick green) with convergence to time-averaged limiting (uniform in this case) 
distribution (thick red) for a seven-cycle with symmetric initial coin state and Hadamard coin operator. The difference between 
the time-averaged distribution and the limiting distribution is scaled by a factor of 0.25 to make the range of fluctuations 
roughly comparable with the entropy, as shown by the instantaneous entropy and distribution difference (thin lines). 



evolution will be more complicated than a line. Cycles with odd or even values of N give different results. For an 
even cycle, only even (odd) positions are occupied after an even (odd) number of time steps, but for an odd cycle, 
after the first {N + l)/2 steps, both even and odd positions are occupied at the same time. We use the same coin 
operator as for the walk on the line, given by equation ([p. We observed that the entanglement of the coin, apart 
from the particular cases identified by Tregenna et al [20| in which the walk is periodic, shows no regular pattern, 
being apparently chaotic. We give two examples illustrating this in figure 1121 Since the entanglement follows the 
instantaneous state of the system, it does not tell anything useful about the mixing properties of the time-averaged 
distribution. We also calculated the time-averaged entanglement (shown in figure [T^ and this appears to converge to 
a steady value at roughly the same rate as the distribution converges to its limiting distribution. This is shown for a 
seven-cycle in figure fT3l 
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FIG. 14: A six-cycle where the numbers ±1, and ±3 label the directions which the particle can take. This is an example of a 
complete bipartite graph, as redrawn in equivalent form on the right. Each blue vertex is connected to every one of the red 
vertices and vice versa; all the edges cross between the two sets of vertices (divided by the dashed line). 
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FIG. 15: (Left) Entanglement Ec{t) for the graph of figure [141 using a DFT coin operator and a symmetric initial coin state. 
(Right) Entanglement Ec{t) for 10-cycle with diagonals and a DFT coin operator. 



B. Cycles with diagonals 

Next we consider the case where there are three or more possible directions that the particle can take. As a 
generalisation of the cycles, we take the case (for N even) where opposite vertices of the A^-cycle are joined to give 
three possible directions for the particle. Figure [Ml explains how this works. The third possible path of the particle 
connects the original position {x = 0) to the opposite position {x — 3) of the cycle, and likewise for the other two 
opposite pairs. The edges of the cycle need to be consistently labeled (see 0, H^])) and here we have chosen to label 
each end of the edges with either —1, 1, —3 or 3, such that adding the vertex and edge label gives the vertex label at 
the other end of that edge. This means that as the particle traverses the edge, the sign of the coin state must flip, so 
we adjust the conditional shift operation to act as 

S\c,v) ^ \- c,v + c{inod N)), (16) 

where v e {0 . . . {N — 1)} is the vertex and c is the coin state, compare equation Although this means we are 
using more coin states than the degree of the graph (four instead of three) , at any single vertex only three of the four 
coin states are actually used, and we pad the coin operator with zeros (one on the diagonal) for the unused coin state 
so it operates correctly (as a Grover or DFT coin) on the three-dimensional subspace. More details on how to do this 
for the general case of a graph with vertices of various degrees can be found in [23| . 

The evolution of the walk is determined by the coin operator: we used the d = 3 Grover and DFT coins, equations 
(fTO|) and ([T2|) respectively. For = 6, the entanglement is periodic with period four, as shown in figure[T5l (Left). For 
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all other values of N that we tested, the TV-cycles plus diagonals showed no regularity, and the entanglement followed 
a complicated pattern, as illustrated in figure [12] (Right) for the case of = 10. 



C. Complete bipartite graphs 



The six-cycle with opposite vertices linked is an example of a complete bipartite graph. In a bipartite graph the 
vertices can be divided into two distinct sets such that every edge connects between the two sets. A bipartite graph 
is complete if, in addition, each vertex in one set is connected to every vertex in the other set. These conditions fix 
the relationship between the number of vertices and the degree of the graph: a complete bipartite graph of degree 
d has N = 2d vertices, we denote it by Kd^d- These graphs can be obtained from the A^-cycle by adding edges 
between every pair of vertices (i, j) in which i is odd and j is even. We can label t he directions from each vertex as 
±1,±3, . . . ,±((i — 1), as shown in figure [T4l for d = 3, i.e., N = 6. lAhmadi et al.l 30l | studied the continuous time 
quantum walk on similar graphs, looking for instantaneous mixing (i.e., a uniform distribution obtained at a particular 
instant in the time evolution of the quantum walk). They found only a small number of examples of instantaneous 
mixing, for regular complete and cyclic graphs with no more than four vertices. This is in sharp contrast to classical 
random walks, which approach a uniform distribution as they evolve on all well-behaved graphs. Periodic behaviour 
is also a property of the instantaneous distributions, though a slightly less stringent requirement than instantaneous 
mixing. 

We can show analytically that the Grover walk on K^^d is periodic with period 4, by looking at the evolution 
operator [/ = 5 • (1 (g) C^^^) in Fourier space. With our chosen edge labeling the shift operator acts as in equation 



16)1 . Taking a FT on the vertex space only. 



1 ^ 



fc=i 



shows that the shift acts as 



1 ^ 



and then 

meaning that S is block-diagonal in the FT basis. Represent it by 
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in which lo — ^^'■1^ ^ and the last entry is in the diagonal for the case where d is odd, otherwise the last block is 
similar to the others. In this basis, the (A^d-dimensional) evolution operator IJ = S ■ C'd*^^) factorises into a 
set of N matrices C/fc each of dimension d, with [7^ = ■ C'^d^. From this one can check explicitly that = 1, 
which means that the walk has a period of 4. Alternatively, by looking at the eigenvalues of Uk we can also verify 



3,...,(<i-l) 



N 



(method used in Q), which 



the periodicity. Uk has eigenvalues ±1 and e with cos 9k = 
reduce to ±i when d = N/2. 

Using numerical simulations with a DFT coin on complete bipartite graphs we found no examples of periodicity. 



WALKS ON THE "GLUED TREES" GRAPH 



We now turn out attention to the special graph used by IChilds et al.l [T| for their algorithm with an exponential 
speed up. An example of this "glued trees" graph with tree depth TV = 4 is shown in figure [TCI fright). At the centre, 
each leaf node has two edges joining it to the leaf nodes of the other tree so, except for the entrance and exit, exactly 
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, Exit 



Entrance 



1 23 4 columns 5 6 789 




Exit 



FIG. 16: Left: regularly glued binary tree graph. Right: "glued trees" graph used in the algorithm of l]. Example shown is for 
A'^ = 4, with 2N + 2 = 10 columns and 2(2''^^^-' — 1) = 62 nodes. The gap between columns 4 and 5 is for clarity in the figure 
and is not significant in the algorithm. 



three edges meet at each node. The problem is to travel via the edges from node to node as quickly as possible starting 
at the entrance and finishing at the exit. The time taken to reach the exit is an example of a "hitting time" (see 
Kempe (sij for definitions and an earlier example of a "hitting time" quantum walk problem). A similar graph with 
a regular join at the leaves is also shown in figure [16] (left). This graph is easy to traverse with a classical algorithm 
because it is easy to identify the middle (in this case because the middle nodes have only two edges joining them, 
but a regular pattern of two edges per node is also classically tractable). The randomly glued edges joining the two 
halves of the graph in figure [TBI (right) disguise the join, and a classical algorithm will get lost at this point, taking 
o n average exp onentially longer to emerge at the exit. 

IChilds et al.l use a continuous time quantum walk for their algorithm. The adjacency matrix A of a graph is 
an N X N matrix with entries Aij = 1 iff there is an edge joining nodes i and j, all other entries in A are zero. For 
an undirected graph (edges can be traversed in either direction, from i —> j and j — > i), A is symmetric. Thus it can 
be used to form the Hamiltonian for the quantum walk: 

*|(^IV'(t)>=E(^|H|y)(y|V(i)), (21) 
y 

with H = 7A where 7 is the transition probability, and where x, y are nodes in the graph. The solution may be 
written 

|^(t))^e-^^*|V'(0)), (22) 

though of course actually calculating it for specific instances of A and |V'(0)) is in general a nontrivial task. The 
proof [l| that the quantum walk is exponentially faster than any classical algorithm involves detailed consideration of 
oracles, colourings and simulation of a continuous time quantum walk on a discrete gate-model quantum computer. 
We will not need to discuss these details here. Figure [17] shows an example of the propagation of the continuous time 
walk through the glued trees graph in terms of the column positions shown in figure 1161 

A discrete time walk can also traverse this graph efficiently if a three dimensional Grover coin is used H ^, . An 
example of the propagation using a Grover coin operator is shown in figure [T8l As noted bv lTregenna et al.. [20], the 
fast hitting time obtained with a quantum walk is highly sensitive to the symmetry of the problem: for quantum 
walks starting at a node other than the entrance, the exit becomes exponentially harder to find and the quantum 
walk does no better than a classical algorithm. 
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FIG. 18; Propagation of the discrete quantum walk using the Grover coin operator on a glued trees graph for tree depth 
= 10, run for 44 time steps. 
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Entrance 




FIG. 19: Detail of a graph with higher branching rate at each node. The branching rate in this example is B — 10, and the 
degree of each node (except the entrance and exit nodes) includes the edge from the parent node soisd = _B + l = ll. 

IXregenna et al.l also noted that that if a DFT coin operator is used instead of a Grover coin operator on the "glued 
trees" graph, the quantum walk stays near the starting point and does not spread out even as far as a classical random 
walk. We consider how this picture changes if we increase the branching rate of the trees that form the graph. In 
other words, we make a similar graph using a pair of trinary trees (branching rate 3) and quaternary trees (branching 
rate 4) and so on for arbitrary branching rate B. This is illustrated in figure fT9l 

A. Mapping to a walk on the line 

Despite the random connections in the centre, the "glued trees" graph with any branching rate is still highly 
symmetric. Provided the initial state used at the entrance node respects the symmetry, the whole quantum walk 
process can be mapped to a walk on a line corresponding to the column positions shown in figure 1161 with different 
biases in the probabilities for moving right or left at each step. The mapping for the continuous time walk is given in 
[l[ for branching rate B — 2. For arbitrary branching rate this generalises to give a Hamiltonian for column positions 
j with non-zero matrix elements 

OTO + l) = {^\°!i.<^'^<^'^^^ (23) 

and (j + l|^f|j) = {j\H\j + 1). We have a choice for the hopping rate 7. In order to make a fair comparison between 
different branching rates B, we take 7 — B~^^^. This makes the non-zero matrix elements of H unity, except at 
the glue, corresponding to unit hopping rate between column positions in the mapped-to-line version of the quantum 
walk, for all choices of N and B. 

For the discrete time quantum walk the mapping is coin specific, and only works when the coin operator preserves 
the symmetry of the graph so the amplitude of the quantum walk is the same on all nodes of each column. To 
perform the mapping for the Grover coin operator given by equation ([9]) , we consider one step of the evolution at a 
typical node in the left hand tree. Our notation is shown in figure [201 We assume that Ci = C2 • • • = Cb = C and 
Di = Z?2 • • • = Db = D. Applying the Grover coin operator gives two relationships between the incoming and the 
outgoing amplitudes, 

D ^ -^P+^C 
B+l B+1 

We then require B^ = B?j and B^ = Li] for the incoming probabilities on the full tree and on the line, 

similarly B^Q^ = Lq and J2a = outgoing probabilities. This gives the relationship between the 
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B-^ nodes in column j 




FIG. 20: Typical node in left hand tree before and after the coin operator is applied, showing our notation for mapping the 
walk to a line, P, Q, {Ci}, {Di} Ri, Ro, Li, Lo are amplitudes. 

amplitudes on the full tree and on the line (independent of choice of coin operator) , 

P = B~^/^Rj C - B-^^+^^/^Li 

Q = B-^/^Ro D = B-^^+^^'^Lo (25) 

Substituting into equation (|24p gives the 2x2 operator for a quantum walk on a line corresponding to the Grover 
coin operator, 

Ro 
Lo 

The right hand tree is a mirror image of the left hand tree, so for that we just exchange L and R in the above 
equations. Unlike the continuous time walk, there is nothing special for the random edges in the glue, the amplitude 
for traversing each edge is determined by the coin operator at the nodes, and this coin operator is different for each 
half of the tree. 

For the entrance and exit nodes we need to consider the most general choice in d dimensions that respects the 
symmetry of the graph. As explained in [2^, although the roots of the tree are only of degree d — 1 = i?, we can 
pad the extra coin dimension with a piece of the identity operator so we have d dimensional coins throughout the 
walk. This has the added effect of allowing an extra arbitrary phase, the most general form of the coin operator at 
the entrance and exit nodes takes the form 

Cend = ^ Q -[^ ^ , (27) 

where C^^^ is a suitably symmetric coin operator of dimension B = d — 1, with an extra arbitrary phase <j), and the 
boldface zeros fill in the row and column to make Cend have dimension d. We then map Cend to a walk on the line, 
which gives a reflection with a phase shift (that may be chosen independently at the entrance and the exit), 

Ro = e^'^^Li Entrance 

Lo = e^-^^i?/ Exit. (28) 

This mapped-to-the-line operator already tells us how the Grover coin operator will behave in the limit of large 
branching rate. Equation becomes —idy as B — > oo, for which the walk on a line simply oscillates between the 
initial and neighbouring nodes, making no further progress along the line. Thus we expect the probability of reaching 
the exit to fall towards zero as B ^ oo. 

We can relate the entanglement between the coin and the position on the line to the entanglement between the 
coin and the position in the walk on the full graph. To calculate the entanglement, we calculate the entropy of the 
reduced density matrix for the coin (by tracing over the position) and then obtain the eigenvalues of this density 



2^/B 
B + 1 
1 - B 
B + 1 



B - 1 
B + \ 



Ri 



(26) 
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FIG. 21: Variation of exit probability with tree depth on the "glued trees" graph with branching rate B = 2 for the continuous 
time quantum walk (black) and the discrete quantum walk with a Grover coin using two different reflection phases as detailed 
in the key (red and blue). 



matrix, from which the entropy is calculated according to equation ([6]). For a two-dimensional matrix the eigenvalues 
are the solution of a quadratic equation, for a d-dimensional coin the eigenvalues are the solution of a polynomial of 
order d. By calculating the reduced density matrices for the full walk and the mapped-to-line versions, and comparing 
the polynomials we find, for example, for the Grover coin case, the polynomials are related by 

f^^lix) = (-x)^-i(:r2 _ ^ ^ ^(^^^^ _ ^1)) ^ (-x)^-V4n^(^), (29) 

where as, cb, ds are time dependent coefficients. The roots of this polynomial give the eigenvalues and thence the 
entropy of the reduced density matrix, which measures the entanglement. Thus we see that the entanglement is the 
same between the coin and position in the mapped-to-line walk as in the corresponding full walk. We studied the 
entanglement for the walk on the glued trees graph, but it did not provide useful indications of the progress of the 
walk, so we do not present any of these results here. 



B. Variation with tree depth A'^ 



We first restrict our attention to the original problem with branching rate B = 2, and consider what happens when 
the tree depth N is varied. The quantum walker does not arrive the exit with certainty after a short number of time 
steps. However, the probability of finding it at the exit node shows a clear peak after roughly the number of steps it 
takes to walk deterministically from entrance to exit by the shortest route. If this peak is large enough (polynomial in 
N and B), one can simply repeat the quantum walk a polynomial number of times to increase the success probability 
to near certainty (probability amplification). 

We studied both continuous time and discrete time quantum walks using simulations with the dynamics mapped 
onto a walk on the line. For the continuous time walk the results are s traightforwar d, given by the black line in figure 
[5T] This is well- fit for large N by P(exit) ~ iV~^/'^, as predicted bv lChilds et al.l [l|, who derive the scahng of the 
Green's function to be '-^ N~^^^ at the peak. 

For the discrete time walk, as noted in the previous section, we have a choice about what to do at the entrance 
and exit where the nodes are of degree two rather than three. We tested two examples of choices of phase shifts, 
equation ((27)) . using the same phase at both entrance and exit. The 2D version of the Grover coin operator is (Jx, 
which produces no phase shift when mapped to the line. Another possible symmetric 2D coin operator is 

of (;.;). (30, 

Mapped to the line, this produces a phase shift of (1 -I- z)/\/2 = e*'^/^. The different phases give different results due 
to the different interference effects. Results for both are also shown in figure [2T] They produce a slightly higher exit 
probability than the continuous time walk for most values of N, but follow the same scaling of P(exit) ~ 7V^2/3 
(numerical data up to iV = 10000 not shown i n figure!^ was u sed for the fitting). This corresponds with the analytic 
solution for a simple walk on a line found by lAmbainis et al.l [isj , which gives the scaling of the peak amplitude as 
~ N~^/^. The mapped-to-the line quantum walk is not exactly the same as the simple walk on an infinite line, but 
the scaling of the peak should be similar, up to the first reflection at the exit. This scaling, being polynomial in 
N, does not affect the exponential speed up of the algorithm. We can apply probability amplification efficiently in 
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FIG. 22: Variation of the first peak in the exit probability for branching rates up to _B = 50, and tree depths up to A'^ = 50 
using a Grover coin operator. 



polynomial time and still beat the exponentially slow classical algorithm. However, the difference between the two 
different choices of 2D coin operators is instructive, suggesting that in other algorithms an appropriate choice might 
significantly improve the success probability. 



C. Variation with branching rate 



We now consider how the exit p robabilit y varies if the branching rate of the two trees forming the graph is varied. 
For the continuous time walk, Chil ds et al.l [l[ provide most of the calculations, they give the transmission amplitude 
through a "defect" (i.e. the glue) of 

^/ X 2iy/BsmK 

(B-l)cosK+(B + l)sinK' ^ ' 

(in our notation) where < k < tt is the momentum of the quantum walker. This has a broad peak at k = 7r/2, 
where the transmission probability |T(7r/2)p — AB /{B + 1)^. Thus for large B we expect the success probability to 
scale as 1/B, and we observed this numerically. 

For the discrete time walk the results are shown in figure [22l The peak exit probability closely follows that of 
the continuous time walk scaling as 1/B. Thus both continuous time and discrete Grover coined walks still beat a 
classical algorithm, which takes time proportional to the number of nodes, i.e. T ~ B^ , to find the exit with better 
than exponentially small exit probability. Combining the scalings for depth and branching rate, we have these two 
versions of quantum walks finding the exit in time linear in TV with probability proportional to 1/{BN'^/^). 



D. DFT coin operator on glued trees 



We now consider the quantum walk with a DFT coin operator on the glued trees graph. The DFT coin operator 
does not have the right symmetry to map to a walk on the line for randomly applied coin labels, but if we apply 
a regular coin labeling that singles out the "root" direction to have coin label zero the mapping to the line can be 
carried out successfully. This is "cheating" as far as application to traversing the graph is concerned, since such a 
labeling would enable a classical algorithm to deterministically find the exit. Nonetheless, the branching trees that 
make up this graph appear widely in physics and computer science, so other applications that could use this regular 
labeling may arise. Note that the coin labeling in [2Q| was consistent in that each end of each edge had the same 
label, but therefore could not always label the root direction as zero, and in fact made random label choices once the 
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FIG. 23: Propagation of the discrete quantum walk using the DFT coin operator on a glued trees graph for tree depth A*" = 10 
and branching rate B = 10, run for 44 time steps. 



consistency condition was satisfied. Hence they did not find any advantage to using DFT coin over classical. With 
the regular labeling, the mapping works because the phase factors all neatly cancel out, as we now show. We start 
by assuming Ci = C2 ■ ■ ■ — Cb — C, and show that after applying the DFT coin operator as in equation p3|) for 
d = B + 1 we have Di — D2 ■ ■ ■ = Db = D. Thus if we start in a symmetric state, this will be preserved under the 
DFT coin operator. Using the notation in figure \20\ 



Q = 

Di = 
D2 = 

Db = 



1 




- 1 


1 






- 1 


1 






- 1 


1 






- 1 



{P + BC} 



.2B 



)C} 



^^\p+{lu^+uj^^ + ---+uj^')c] 

; + i I J 



(32) 



where w = e^Tri/d ^j^g complex d'th root of unity. Each of the sums of powers of uj in the expressions for Di to Db 
evaluates to —1, hence we have Di = D2 ■ ■ ■ = Db = D as claimed. Substituting for P, Q, C and D from equation 
then gives 



Ro — 
Lo = 



B + 1 
1 



B 1 



VB + T 



Ri 



B 



B + 1 



-L, 



(33) 



for the 2x2 operator equivalent to the DFT coin operator. For B ^ 00, equation (j33p becomes the identity, which 
corresponds to a walk that steps deterministically along the line, reaching the exit in the shortest possible time. For 
intermediate values of B, the DFT coined walk partly steps across the graph, and partly remains close to the entrance: 
an example is shown in figure [23] for a "glued trees" graph with branching rate i? = 10 and depth iV = 10. 

Again we look at the first peak in the exit probability, which occurs after time ~ (27V + 2). The variation of this first 
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FIG. 24: Variation of the first peak in the exit probability for branching rates up to -B = 50, and tree depths up to A'' = 50 
using a DFT coin operator. 

peak in the probability with branching rate and tree depth for N,B < 50 is plotted in figure [24l For B — N, there is 
a wide peak; on the B > N side an initial peak rises towards unity, while on the B < N side the probability tails off 
in a series of smaller arches. For B — 2, the quantum walk with this coin operator has a noticeably lower probability 
of finding the exit after linear time than the Grover coin or the continuous time quantum walks. Nevertheless, there 
is still a peak in the probability of finding the walker at the exit, and our numerical results show this peak scales as 
N~^^^, the same as the Grover coined walk and the continuous time walk, but with a constant prefactor ~ 30 times 
smaller. 

We can understand this prefactor qualitatively as follows. The mapped-to-the-line version of the coined quantum 
walk is essentially a generalised Hadamard coin walk on finite line segments f22|. Such walks on the infinite line all 
have the typical double-peaked shape, with the ratio of the left and right peak heights determined by the coin biases 
and the initial state. The mapped-to-the-line versions of the Grover and regular DFT coins have the minus sign in 
opposite positions in the off-diagonal, equations (I26|) and (I33p . Combined with the rightward moving initial state, for 
the first half of the walk the largest peak is thus on the leading(trailing) edge for the Grover(DFT) coin operator. 
Moreover, the trailing edge DFT peak is larger than the leading edge Grover peak. For the second half of the walk 
the coin operators are transposed, and the traihng(leading) edge is largest for the Grover(DFT) coin operator, but 
this is combined with different "initial" states as the quantum walk crosses the glue from the first half of the walk, 
so the two effects don't just cancel out. The cross-over point beyond which the regular DFT coin performs best is 
roughly S > 10 for any tree depth. 

VI. SUMMARY AND DISCUSSION 

Even for the simple case of a coined quantum walk on the line we find that the entanglement between the coin and 
the particle position shows surprisingly rich behaviour. The entanglement oscillates around an asymptotic value that 
is determined by the bias in the coin operator, and the rate of convergence depends the symmetry in the distribution 
of the quantum walk, which in turn is determined by the initial state of the coin (assuming the quantum walk starts at 
the origin). More symmetric distributions (about the origin) converge faster than asymmetric distributions in which 
there is a bias favouring finding the particle in the positive or the negative half line. 

For lattices in two spatial dimensions, similar rich patterns of behaviour are observed, the rate of convergence of the 
entanglement correlates with the extent to which the quantum walk spreads out over the grid, with fast convergence 
indicating good spreading behaviour, and slow convergence for distributions that bunch around one point, such as 
the origin. We find that for lattices in two spatial dimensions but of higher degree c? > 4, a DFT coin operator is 
better at spreading over the graph than a Grover coin operator. This can be explained by the bias in the Grover coin 
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operator: above degree d = 4, the Grover coin operator favours returning along the path it arrived from rather than 
picking a new direction, thus tending to keep the walk near where it started. The DFT coin operator is unbiased for 
any degree, so the spreading depends only on the relative phases, which can be controlled by the coin labeling. 



For finite graphs, there are a few periodic examples, corresponding to those reported in [20], otherwise the entan- 
glement is random even for very simple cases, like the seven-cycle. A quantum walk on a complete bipartite graphs 
using a Grover coin operator is always periodic with period four, and we have shown this analytically. For the DFT 
coin, the quantum walk on complete bipartite graphs was not periodic for any cases we tested numerically. 

For the glued trees graph and its generalisation to higher branching rates, the first peak in the probability for finding 
the quantum walker at the exit node occurs in linear time and scales with tree depth as N~^^^ for both Grover and 
DFT coin operators, matching that found for the continuous time walk bv lChilds et al.l This is in contrast with 
classical algorithms, which require exponential time to find the exit with better than exponentially small probability. 
However, for increased branching rate, while the continuous time walk and the Grover coined discrete time walk both 
scale as B~^, a DFT coin with regular coin labels has a success probability approaching unity ioi B > N >> 1, 
and beats the continuous time and Grover coined walks for i? > 10. A regular coin labeling is not permitted in the 
algorithmic context of but branching trees occur in many other contexts, some of which may permit a regular 
coin labeling such that this DFT coin property could prove useful. 

Another observation from this work is that since the coin needs to have as many dimensions as the maximum degree 
of the graph, it allows a relative phase to be introduced at vertices of less than the maximum degree, for example, at 
the entrance and exit of the "glued trees" graph, see figure 1211 The quantum walk may exhibit different properties 
depending on the choice of relative phase, which provides extra options for tuning the quantum walk. 

Overall, the discrete quantum walk offers many options for tuning the properties of the quantum coin operator to 
fit the problem to be solved. These possibilities are not available in the continuous time quantum walk, unless one 
adds extra degrees of freedom [l§| (whereupon it becomes exactly equivalent). Further work could investigate how to 
optimise the uniforrnity of spreading over these higher dimensional lattices, for example, by applying small amounts 
of decoherence (see [20]) and by adjusting the coin operator: we do not know whether the DFT coin operator is the 
optimal choice for a general rf-dimcnsional graph for > 4. 
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